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ABSTRACT 

From a study of the integrated light properties of 200 globular clusters (GCs) in M31, Strader 
et al. found that the mass-to-light ratios are lower than what is expected from simple stel¬ 
lar population (SSP) models with a ‘canonical’ stellar initial mass function (IMF), with the 
discrepancy being larger at high metallicities. We use dynamical multi-mass models, that in¬ 
clude a prescription for equipartition, to quantify the bias in the inferred dynamical mass 
as the result of the assumption that light follows mass. For a universal IMF and a metallicity 
dependent prevent day mass function we find that the inferred mass from integrated light prop¬ 
erties systematically under estimates the true mass, and that the bias is more important at high 
metallicities, as was found for the M31 GCs. We show that mass segregation and a flattening 
of the mass function have opposing effects of similar magnitude on the mass inferred from 
integrated properties. This makes the mass-to-light ratio as derived from integrated properties 
an inadequate probe of the low-mass end of the stellar mass function. There is, therefore, no 
need for variations in the IMF, nor the need to invoke depletion of low-mass stars, to explain 
the observations. Finally, we find that the retention fraction of stellar-mass black holes (BHs) 
is an equally important parameter in understanding the mass segregation bias. We specula¬ 
tively put forward to idea that kinematical data of GCs can in fact be used to constrain the 
total mass in stellar-mass BHs in GCs. 

Key words: globular clusters: general - open clusters and associations: general. 


1 INTRODUCTION 


The stellar initial mass function (IMF) plays a vital role in astro¬ 
physics and it is an ongoing debate whether the IMF is universal for 


all metallicities 

Bastian et al. 20101: van Dokkum & Conroy 2010; 

Gehaetal, 2013 

). Globular clusters (GCs) are important test-beds 


for studies of the IMF and the present day mass function (MF) be¬ 
cause they are, in most cases, single aged stellar populations with 
a single [Fe/H], and the kinematics of the stars within them is due 
to the gravity of the stars only. Because most of the GC mass is 
in low-mass stars, whereas the luminosity is due to the high-mass 
stars, the GC mass-to-light ratio T y = M/ Ly is consider to be 
a good probe of the shape of the stellar MF. Here Ly is the V- 
band luminosity and M is independently found from a measure¬ 
ment of the cluster’s radius and velocity dispersion. Recent obser¬ 
vational advances have allowed us to measure the structural pa¬ 
rameters and kinematical properties of GCs in external galaxies 
up to distances of several tens of Mpc ISmith & Gallagherll200ll : 

iBastian et alJ2006l : lBrodie & Straden2006l) . spanning a large range 

of environments, thereby potentially increasing the applicability of 
T y as a probe of the IMF. 

For Galactic GCs, a slight positive correlation between Ty 
and M was found by early work by iMandushev et al. | « 199ll) 


and IPrvor & Mevlanl dl993l) . A more significan©correlation was 
found from a larger sample of G Cs by iKimmig et al.l d2014h . 
iMcLaughlin & van der Mareil ( t2005t) found an average (Tv) — 
1.45 for Milky Way GCs, which is slightly lower than the SSP pre¬ 
diction of Tv — 2, appropriate for low metallicities. This small 
difference and the trend with M has been attributed to the pref¬ 
erential escape of low-mass stars over the tidal boundary, reducing 
M, but hardly affecting Ly (e.g. lKruiissei3l2008l) . 

I Strader et al .1 1201 ll hereafter S11) published the T v values 
of 200 GCs in M31, spanning a range of metallicity of —2 < 
[Fe/H] < 0. They found that at low metallicities the average Tv 
agrees well with the values from simple stellar population (SSP) 
models, assuming a Rroupa IMF, albeit with a scatter of a factor 
of two. For increasing metallicities, a growing difference between 
the dynamical Tv and the SSP models was found, with Tv falling 
below the SSP values, by up to a factor of 3 at [Fe/H]_^ 0. Similar 
results were recently obtained for Milky Way GCs ( Kimmig et al J 
120141) . These results are hard to explain by the preferential escape 
of low-mass stars, because it implies that almost all metal-rich clus¬ 
ters have lost a significant fraction (~ 60% — 70%) of their initial 
mass. 

An alternative interpretation is that the IMF gets more bot- 
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tom light (fewer low-mass stars) with increasing [Fe/H]. Before 
accepting the possibility of a varying IMF, we need to establish 
that there are no metallicity dependent biases as the result of dy¬ 
namical evolution of GCs in the Tv measurements. A possible 
bias is the effect of mass segregation on the measurements of the 
GC radius and velocity dispersion. Mos t GCs are older than their 
h alf-m ass relaxation time-scale r r h (e.g. lHeno ni l 196 ll : [Gieles et al] 
1201 lh . This means that they had time to evolve towards equiparti- 
tion and, as a result, the most massive objects (stellar remnants and 
turn-off stars) are more centrally concentrated and move with lower 
velocities. For resolved GCs this effect can be taken into account 
by fit ting multi-mass models (e.g. IPaust et aT]|2() I (jl : ISollima et al.1 
I2012T) . which are distribution function based model s that inclu de an 
approximate prescription fo r equipartition jDa Costa & Freemanl 
1 19761 : iGunn & Griffir3l 19791 hereafter GG79). Flowever, in extra- 
galactic samples where the individual stars are not resolved, this 
is harder to do (although not impossible), and it is therefore often 
assumed that light follows mass. 

For a universal IMF, the stellar mass function at an age of 
12Gyr does depend on the metallicity of the stars because of dif¬ 
ferences in stellar evolution, and as a result we may expect that 
clusters with different metallicities, have different observable prop¬ 
erties. The idea that a metallicity dependent bias in the observation- 
ally derived cluster properties may exist is not new. |jor dliil 120041) 
considered multi-mass King models, with mass bins appropriate for 
different metallicities. Fie showed that as a result of mass segrega¬ 
tion, clusters with the same half-mass radius ru and the same age, 
have smaller half-light radii in projection, rh p ,£, at higher metal¬ 
licities. This is because at high [Fe/H] the turn-off mass is higher 
and as a result these stars and the bright evolved stars, which have 
similar masses, are more centrally concentrated than in metal-poor 
clusters. Jordan proposed that this effect is responsible for the ob¬ 
served size difference between metal-rich and metal-poor GCs (e.g. 
lLarsen et alJhoOlT ). An apparent size difference between metal- 
rich and metal-po or GCs was also found in JV-b ody simulations 
(Sip pel et ai]|2012l) and Monte Carlo simulations jDownin mu; 
of clusters with different metallicities. 

In this Letter we use multi-mass King models to quantify the 
mass segregation bias in Ty as a function of [Fe/H] and for dif¬ 
ferent assumptions for the retention of remnants. In §[2]we model 
stellar mass functions for different [Fe/H] and describe the dynam¬ 
ical models. In §[3]we present the results of a comparison to the data 
by SI 1 and our conclusions and discussion are presented in §[4] 


2 DESCRIPTION OF THE MODELS 
2.1 Stellar mass functions at an age of 12 Gyr 


Here we describe how we compute the stellar mass components 
for clusters with different metallicities. We create N p = 10 4 star 
particles with initial masses m i equally spaced in log mi between 
a lower mass of m\ 0 = 0.2 M@ and an upper mass of m up = 
100 M@. We consider a power-law (‘canonical’) IMF 


d N 
dm; 


= Am,: 


( 1 ) 


where a = 2.35 and A is a constant that normalises the to¬ 
tal number of stars in the cluster: A = (1 — a)/(m \“ Q — 
mj L 0 _ “), f° r a ^ 1. Every star particle is given a weight w = 
dm[ _ “(lnm U p — lnmio)/JV P , such that = 1 and 

vjjrri\j ~ 0.68 Mq is the mean mass of the IMF, which 


is similar to the value for more realistic IMFs that roll over at low 
masses to flatter slopes and are truncated at m\ 0 = 0.1 Mq. For 
our purpose this suffices, because the behaviour of the multi-mass 
models is sensitive only to the total mass in low-mass stars, and in¬ 
sensitive to the mean mass of the low-mass stars (see Section 5 of 
GG79). Tests confirm that our results are insensitive to this choice. 
To study the effect of mass segregation in a system with a bottom 
light IMF, or a cluster that is depleted in low-mass stars because of 
dynamical evolution, we also consider a double power-law (‘flat’) 
IMF 

d N_ 

dm; 

where B = A to make the function continuous. We normalise the 
high-mass end in the same way as for the ‘canonical’ IMF (i.e. the 
same A), so the reduction in mass due to flattening can be found 
straightforwardly. 

We evolve every star particle to an age of 12 Gyr using SSE 
d Hurley et al.ll200ol ) for 11 metallicities, linearly spaced in —2 ^ 
[Fe/H] < 0. We split the data into bins corresponding to the stel¬ 
lar types at 12 Gyr: main sequence stars (MS), evolved stars (EV), 
white dwarfs (WD), neutron stars (NS) and black holes (BH). The 
logarithmic spacing] in im and the fast SSE code allow us to effi¬ 
ciently model MFs of old stellar populations that are well sampled 
at all evolutionary stages and hence present day masses m. 

The resulting MFs for a metal-rich and a metal-poor cluster 
with a ‘canonical’ IMF are shown in Fig.Q] From this figure we see 
that at [Fe/H] = —2 the turn-off mass is slightly lower (0.8 M 0 ) 
than at [Fe/H] = 0 (1 M 0 ) and that the average mass of WDs is 
higher at low [Fe/H] (0.77 M 0 ) than at high [Fe/H] (0.65 M 0 ). 
Also, the contribution of remnants to the total mass is higher at low 
[Fe/H]. Although these differences are small, they are important 
for the distribution of the visible stars, as we will show in the next 
section. 

For our dynamical models we create 13 mass components: we 
split the MS stars into 5 bins that are equally spaced in log(m), the 
WDs into 5 bins that are equally spaced in m, the EV stars, NSs 
and, finally, the BHs. With the [Fe/H] dependent mass functions in 
place, we can now describe the dynamical multi-mass models. 


B, mio < m i ^ 1 M 0 , 

Amf a , 1 M 0 < rrii < m up , 


( 2 ) 


2.2 Multi-mass models 

In this section we describe the dynamical models that we use to 
quantify the effect of mass segregation on the inferred dynami¬ 
cal mass from observations. For a self-gravitating system in virial 
equilibrium, the mass M can be written as a function of the mass- 
weighted mean-square velocity of the stars, ( v 2 ), and the virial ra¬ 
dius of the cluster, r v , as M = 2 (v 2 )r v /G. Here r v is defined as 
r v = —GM 2 /(2U), where G is the gravitational constant and U 
is the total potential energy of the cluster. Our focus is on extra- 
galactic cluster samples for which individual stars are not resolved, 
so we write the virial relations in terms of the luminosity-weighted 
mean-square velocity in projection {v p )l and the radius containing 
half the total luminosity in projection rh p ,L (or, ‘effective radius’) 


1 Approximating the mean mass in a bin by the central logarithmic value 
is only exact for a power-law with index —1.5. From a convergence test we 
find that for 10 4 stars our model properties such as the total M of the MF 
are accurate to about 10“ 5 . 
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Figure 1. Mass functions for two of our evolved stellar populations: one 
metal poor (top panel), one metal rich (bottom panel). Each distribution 
is split into the five stellar types: main sequence stars (MS), evolved stars 
(EV), white dwarfs (WD), neutron stars (NS), and black holes (BH). The 
numbers are the fraction of the total mass in each component. 


M - 2 71 

rj r rj v G 


(3) 


Here we introduce three factors 77 , r] r and ri v to convert 3- 
dimensional properties to observable quantities, namely 


( v 2 ) fh r v 
{Vp) ?'h P r h ’ 


r/r = 


^hp,L 

rhp 



(4) 


Here { Vp ) is the mass-weighted mean-square velocity in projection 
and rhp is the radius containing half the mass i n projection. The 
value of r/ depends on the model, and for a lPlumme ^ fi9nh model 
with an isotropic velocity distribution p ~ 3 x (4/3) x (5/4) = 5, 
and this value is relatively insensitive to t he choice of the model 
(Spitzer 1987; Portegies Zwart et al. 2010). From now on we as¬ 
sume that uncertainties in 77 are small and can be ignored. 

The factors r/ r and T] v equal unity if light follows mass, i.e. 
if Tv is the same everywhere in the cluster. If this assumption is 
made for clusters for which this does not hold, then the product 
r/ r ri v is the error that is made in the mass estimate. iHurlevl 12007}) 
finds from TV-body simulations that r/ r ~ 0.5 which implies an 
error in M of a factor of two, confirming that mass segregation is 
an important effect to consider. 

We derive the bias factor r) r riv from isotropic multi-mass King 
models as described in detail in GG79. The distribution function of 
each component rrij (section ETt is 


fj(E)=A j 


exp 




(5) 


Here E = | v 2 + </>, and tj> is the (lowered) specific potential and cry 
is a scale velocitjf] that is chosen to depend on my as cry oc 

From the lMarastoni d2005t) SSP models we find that about 65% 
of the U-band light comes from EV stars and 35% from the MS 


2 We note that only in the limit of infinite concentration <ry equals the cen¬ 
tral ID velocity dispersion r/[/ s of bin j. For realistic concentration values 
the ratio r’/j lH /rj y < 1 and is mass dependent. Multi-models models are 
therefore not in equipartition. In forthcoming studies Isollima et al J|20 l5l : 
Peuten et al., in prep) we show that multi-mass models describe the phase 
space distribution of stars in IV-body systems very well. 




Figure 2. Density and velocity dispersion profiles for two of our GC models 
with Wo = 7. Although our multi-mass models have 13 bins, for clarity we 
have only plotted the middle mass bin values for MS and WD. Solid lines: 
[Fe/H] = —2. Dashed lines: [Fe/H] = 0. 


stars. For a typical mass-luminosity relation for MS stars we esti¬ 
mate that about 70% of the MS luminosity is emitted by stars in the 
mass range of EV stars (> 0.65 M 0 ). We therefore assume that the 
dynamical properties of EV stars are the ones that one infers from 
integrated light studies. 

In Fig. [2] we show the resulting density and velocity profiles 
for a metal-rich and a metal-poor cluster without BHs. For both 
clusters the EV stars are more centrally concentrated and move with 
lower velocities than the MS stars which implies that both r/ r < 1 
and T] v < 1. The underestimation of M is more severe for the 
metal-rich cluster than for the metal-poor cluster. This is because at 
high [Fe/H], the turn-off mass is higher, and the remnants are less 
massive (see Fig.QJ, allowing the bright stars to segregate further to 
the centre of the cluster (see also Jordan 2004; Sippel et all 20121) . 

We now need to make a choice for the dimensionless central 
potential Wo- We choose Wo such that the EV stars have a ra¬ 
tio of the core radius ro (or King radius) over rh that is equal to 
what is found for the M31 GCs by SI 1. We obtained the values 
for ro/rh p ,i of the GCs in M31 from Sll and found that there is 
no correlation between ro/rh p ,z, and [Fe/H]. We therefore adopt 
3 values: the mean observed value and the values corresponding to 
the mean and the plus and minus la standard deviation. Because 
r 0 for the individual components has a different meaning then r 0 in 
single mass models because, we decide to define ro as the distance 
to the centre where the density of EV stars has dropped to 1/3 of the 
central density, as is the case in single mass models with Wo > 5. 
Finally, we use rh p ,z, = 0.75rJ/ v , where r|) v is the (3D) half-mass 
radius of the EV stars. 

Using these three ro/rhp.i, values, we produce ten models 
with Wo ranging from 6 to 24, for each [Fe/H]. We then inter¬ 
polate to obtain the values of Wo to obtain the desired ro/fhp.i, 
value. A sample of the Wo values used can be found in Table 1. 


3 RESULTS 

For each [Fe/H], ro/rhp.r and different remnant fractions, we cal¬ 
culate the ratio of the mass as it would be derived from the evolved 
stars, M 0 bs, over the true mass M, i.e. r} r ri v (equation [4}. In this 
calculation, we used exclusively the contribution from the evolved 
stars to derive M 0 b s , in order to be consistent with the observations. 
For the models with a ‘canonical’ IMF M 0 b s /M = M 0 b s /Mssp, 
and for the models with a ‘flat’ IMF we multiplied A/ 0 bs/M 
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Table 1. A sample of the Wq values used for different [Fe/H] and the 
different assumptions for retention of NSs and BHs for a ‘canonical’ IMF. 


[Fe/H] 

r o/ r hp,L 

All BH & NS 

NoBH 

No BH & NS 


0.45 

15.9 

8.86 

8.41 

-2 

0.29 

17.5 

11.1 

10.3 


0.13 

19.8 

14.2 

12.9 


0.45 

23.0 

10.8 

10.5 

0 

0.29 

24.4 

13.4 

12.4 


0.13 

25.8 

16.6 

15.1 


by M/Mssp ~ 0.41, slightly dependent on [Fe/H], to obtain 
Mohs/Mssp. The results are shown in Fig. [3] The three panels 
show the results of three assumptions for the remnant retention 
fraction, from left to right: all BHs and all NSs retained, no BHs 
retained and all NSs retained, no BHs and no NSs retained. The 
results for a model in which all BHs are retained and no NSs are 
retained were found to be almost identical to the model in which 
all NSs and all BHs were retained. 

We divided the inferred Tv of each GC in M31 by its cor¬ 
responding value of the SSP model (Fig. 1 of Sll), based on a 
Kroupa IMF, to obtain M 0 b a /Mssp for M31 GCs, where Mssp is 
the present day mass of a stellar population with a ‘canonical’ IMF. 
We note that there is a slight inconsistency in the comparison be¬ 
tween our models and the data because the Ty values of the SSP 
models are based on a single value for the retention fraction of NSs 
and BHs. However, as we show in Fig.Q]the total mass in NSs and 
BHs never exceeds 9% of the present day M, so this effect is small. 

If all NSs and BHs are retained (left panel), there is hardly 
any bias because of mass segregation and the correct M is in¬ 
ferred from the EV stars. Without NSs and BHs (right panel), 
0.25 < AToba/ATssp ;$ 0.5 for both MFs. The mass segrega¬ 
tion term (r/rVv) is higher when the mass function is flatter, while 
M/Mssp is smaller. The two effects appear to roughly cancel, such 
that the resulting ratio M a b s /Mssp is similar. This means that a 
low M 0 bs/M ssp can not be interpreted as a flattening of the MF, 
nor the IMF, because a mass segregated model with a ‘canonical’ 
IMF gives similar results. 

The two extremes of the retention fraction models (i.e. all the 
BH and NS retained, and none of the BH and NS retained) encom¬ 
pass the majority of the GC values, and it is possible to interpret 
the spread in the data as a spread in BH retention. 

The [Fe/H] dependence in the middle and right panels can 
be explained with some simple stellar evolution arguments. From 
Fig. [I] we see that at lower metallicities the turn-off mass is lower 
and there is more mass in remnants, and the average mass of white 
dwarf mass and black holes is higher (see also the discussion in 
ISippel et al.l2012l) . As a result, in a low-metallicity cluster, the rem¬ 
nants are more efficient in pushing the evolved stars out, which in¬ 
creases their velocities and their half-light radius, such that r/ v and 
r/ r are closer to unity. 

These results suggest that the discrepancy between the ob¬ 
served GC masses of Sll and the SSP models, can be explained 
by metallicity dependent mass segregation effects and a spread in 
the retention fraction of NSs and BHs. Our results suggest that the 
low-mass end of the MF can not be constrained from integrated 
light properties. 


4 CONCLUSIONS AND DISCUSSION 

In this study we show that the discrepancy between the observed 
mass-to-light ratio Tv and simple stellar population (SSP) mod¬ 
els as found in a sample of 200 globular clusters (GCs) in M31 
by Sll, can be explained by mass segregation. We modelled GCs 
with a range of metallicities and with different remnant retention 
fractions. Our key results are: 


• For clusters without stellar-mass black holes (BHs) and neu¬ 
tron stars, the dynamical mass derived from evolved stars ( M a bs ) of 
a mass segregated cluster underestimates the true mass M, and this 
bias is stronger at high [Fe/H] (factor of ~ 4) than at low [Fe/H] 
(factor of ~ 2). This is due to the higher turn-off mass and the lower 
white dwarf masses at high [Fe/H] which cause the evolved stars 
to be more centrally concentrated and move with lower velocities, 
compared to turn-off stars i n low [Fe/H] clusters with the same 
mass density profile (see also]jordan 2004: lSk>Del et alJf2012h . 

• The presence of BHs has a large effect on the derived M 0 bs, 
in that its value is closer to the real M for higher retention frac¬ 
tions. Despite the fact that a large fraction of BHs probably gets 
ejected from GCs as the result of supernova kicks, it is important 
to consider the BHs in dynamical mass modelling of GCs. A small 
remaining BH populatio n can s urvive f or as long as 10 half-mass 
relaxation times (r r h, iBreen & Heggiej|2013l) . The recent discov¬ 
ery of two BH candidates in M22 (Strader et al ■ 20121) confirms the 
need to consider BHs when deriving Ty. 

• For clusters with a ‘flat’ mass function, we find that the effect 
of mass segregation on M 0 b s /M is less important, such that the ra¬ 
tio M l:i bs/M ssp, with Mssp is the present day mass for a canonical 
IMF, is similar to M/Mssp models with a ‘canonical’ IMF. This 
makes Ty insensitive to the low-mass end of the MF. 


Sll find that the discrepancy between Ty and the SSP mod¬ 
els is more important for clusters with lower mass and for clus¬ 
ters with shorter r r h, which was also found for G Cs in the Milky 
Way l lMandush ev et al.l l991 ; Kimmi g et al .120141) and Centaurus A 
(NGC 5128. lReikuba etal] 20071) . This is consistent with the sce¬ 
nario put forward in this Letter, because mass segregation is more 
important for cluste rs with shorter r r h, i.e. those with lower mass. 

[Zaritskv et al.1 l l2014l) report a bi-modality in Ty, after all 
clusters have been aged to a common age of 10 Gyr (Ty,io). They 
determine dynamical masses from integrated light properties of a 
sample of star clusters with different ages and metallicities and in 
different galaxies. The authors interpret their result as evidence for 
two distinct IMFs. Although we do not discuss the mass segrega¬ 
tion bias in their sample in detail here, we caution to invoke the 
need of a varying IMF in a study based on integrated light mea¬ 
surement, before biases as a result of mass segregation have been 
fully quantified. 

The reduced Ty of MW GCs l lMcLaughlin & van der Marell 
120051 ) and other galaxies has been explained by the depletion 
of lo w-mass stars as the result of dynamical ejections l lKruiisseij 
l2008h . This effect is also more important for clusters with shorter 
T r h/lower mass, therefore degenerate with the mass segregation 
bias reported here. We note that preferential depletion of low-mass 
stars requires mass segregation, but not vice versa and it is therefore 
necessary to consider the biases due to mass segregation together 
with the effect of the depletion of low-mass stars. Here we showed 
that the two effects have opposing effects on M 0 b a /Mssp, such 
that Ty is insensitive to the slope of the MF. 

The Ty values of young (few 100 Myrs) massive (> 10® M 0 ) 
clusters in merger remnants are consistent with the predictions of 
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Figure 3. Ratio of inferred mass M 0 b s over SSP mass Msspversus [Fe/H] for a ‘canonical’ IMF (solid lines) and flat IMF (dashed lines). The blue points 
are the corresponding values for the M31 GCs taken from SI 1. The three lines in each of the plots show the results for models designed to Took’ as similar to 
the S11 GCs as possible, using the mean rg/r^p^ value and two a values. For the ‘canonical’ IMFM/Mgsp = 1. and for the‘flat’IMF M/Mssp — 0.41 


SSP models with a ‘canonical’ IMF iBastian et al.l2006t ). Although 
these values were derived from integrated light properties, these 
clusters are significantly younger than their two-body relaxation 
time scale, such that mass segregation is likely not important yet. 

Finally, the strong dependence of M 0 \, a /M on the retention 
fraction of stellar mass BH opens the exciting opportunity to use 
kinematic data of GCs to determine the total mass in stellar mass 
BFls, thereby placing complementary constraints on BH retention 
in Milky Way GCs. 


ACKNOWLEDGEMENTS 

This work was carried out as part of the ISIMA programme on 
Gravitational Dynamics held in July 2014 at CITA, Toronto. We 
acknowledge Pascale Garaud for the organisation, the financial sup¬ 
port and for providing a stimulating environment. MG thanks the 
Royal Society for financial support in the form of a University 
Research Fellowship (URF) and the European Research Council 
(ERC-StG-335936, CLUSTERS). RLS acknowledges the support 
of the Science and Technology Facilities Council (STFC) via the 
award of an STFC Studentship. RLS thanks Anna Lisa Varri for 
her support and guidance throughout this work. We thank Mik- 
los Peuten, Antonio Sollima, Anna Lisa Varri and Alice Zocchi 
for discussions during the ‘Gaia Challenge’ meeting and Diederik 
Kruijssen,Vincent Henault-Brunet, Jay Strader and Nate Bastian 
for comments on the manuscript. Finally, we thank the referee for 
constructive comments that helped to improve the paper. 


REFERENCES 

Bastian N., Covey K. R., Meyer M. R., 2010, ARA&A, 48, 339 
Bastian N., Saglia R. P, Goudfrooij P, Kissler-Patig M., Maraston 
C., SchweizerF., Zoccali M., 2006, A&A, 448, 881 
Breen P. G., Heggie D. C., 2013, MNRAS, 432, 2779 
Brodie J. P„ Strader J„ 2006, ARA&A, 44, 193 
Da Costa G. S., Freeman K. C., 1976, ApJ, 206, 128 
Downing J. M. B„ 2012, MNRAS, 425, 2234 


Geha M., Brown T. M., Tumlinson J., Kalirai J. S., Simon J. D., 
Kirby E. N., VandenBerg D. A., Munoz R. R., Avila R. J., 
Guhathakurta P.. Ferguson H. C., 2013, ApJ. 771. 29 
Gieles M„ Heggie D. C., Zhao H„ 2011. MNRAS, 413, 2509 
Gunn J. E„ Griffin R. F„ 1979, AJ, 84, 752 (GG79) 

Henon M., 1961, Annales d’Astrophysique, 24, 369 
Hurley J. R„ 2007, MNRAS, 379, 93 
Hurley J. R„ Pols O. R„ Tout C. A., 2000, MNRAS, 315, 543 
Jordan A., 2004, ApJ, 613, LI 17 

Kimmig B.. Seth A.. Ivans 1.1., Strader J., Caldwell N., Anderton 
T„ Gregersen D„ 2014, arXiv: 1411.1763 
Kruijssen J. M. D„ 2008, A&A, 486, L21 

Larsen S. S., Brodie J. P., Huchra J. P, Forbes D. A., Grillmair 
C. J., 2001, AJ, 121, 2974 

Mandushev G., Staneva A., Spasova N., 1991, A&A, 252, 94 
Maraston C., 2005, MNRAS, 362, 799 
McLaughlin D. E., van der Marel R. P, 2005, ApJS, 161, 304 
Paust N. E. Q., Reid I. N., Piotto G., Aparicio A., Anderson J., 
Sarajedini A., Bedin L. R., Chaboyer B., Dotter A., Hempel M., 
Majewski S., Marln-Franch A., Milone A., Rosenberg A., Siegel 
M„ 2010, AJ, 139, 476 
Plummer H. C.. 1911, MNRAS, 71, 460 

Portegies Zwart S. F., McMillan S. L. W., Gieles M., 2010, 
ARA&A. 48, 431 

Pryor C., Meylan G., 1993, in Djorgovski S. G., Meylan G., eds, 
Structure and Dynamics of Globular Clusters Vol. 50 of Astro¬ 
nomical Society of the Pacific Conference Series, Velocity Dis¬ 
persions for Galactic Globular Clusters, p. 357 
Rejkuba M., Dubath P, Minniti D., Meylan G., 2007, A&A, 469, 
147 

Sippel A. C., Hurley J. R„ Madrid J. P, Harris W. E„ 2012, MN¬ 
RAS, 427, 167 

Smith L. J., Gallagher J. S., 2001, MNRAS, 326, 1027 
Sollima A., Baumgardt H., Zocchi A., E. B., Gieles M., Henault- 
Brunet V., Varri A., 2015, MNRAS, submitted 
Sollima A., Bellazzini M„ Lee J.-W., 2012, ApJ, 755, 156 
Spitzer L., 1987, Dynamical evolution of globular clusters. 

Princeton, NJ. Princeton University Press, 1987, 191 p. 

Strader J., Caldwell N„ Seth A. C., 2011, AJ. 142, 8 (SI 1) 
Strader J., Chomiuk L., Maccarone T. J., Miller-Jones J. C. A., 
Seth A. C., 2012, Nature, 490, 71 


© 0000 RAS. MNRAS 000, 000-000 



















6 Rosemary L. Shanahan & Mark Gieles 

van Dokkum P. G., Conroy C., 2010, Nature, 468, 940 
Zaritsky D., Colucci J. E., Pessev P. M., Bernstein R. A., Chandar 
R„ 2014, arXiv: 1409.3235 


© 0000 RAS, MNRAS 000, 000-000 


